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Abstract 

In this paper, we propose optimal tests for reflective circular symmetry about a fixed 
median direction. The distributions against which optimality is achieved are the so-called 
fc-sine-skewed distributions of Umbach and Jammalamadaka (2009). We first show that 
sequences of k -sine-skewed models are locally and asymptotically normal in the vicinity of 
reflective symmetry. Following the Le Cam methodology, we then construct optimal (in the 
maximin sense) parametric tests for reflective symmetry, which we render semi-parametric 
by a studentization argument. These asymptotically distribution- free tests happen to be 
uniformly optimal (under any reference density) and are moreover of a very simple and 
intuitive form. They furthermore exhibit nice small sample properties, as we show through 
a Monte Carlo simulation study. Our new tests also allow us to re-visit the famous red wood 
ants data set of Jander (1957). We further show that one of the proposed parametric tests can 
as well serve as a test for uniformity against cardioid alternatives; this test coincides with the 
famous circular Raylcigh (1919) test for uniformity which is thus proved to be (also) optimal 
against cardioid alternatives. Moreover, our choice of fc-sinc-skewed alternatives, which are 
the circular analogues of the classical linear skew-symmetric distributions, permits us a Fisher 
singularity analysis a la Hallin and Ley (2012) with the result that only the prominent sine- 
skewed von Miscs distribution suffers from these inferential drawbacks. Finally, we conclude 
the paper by discussing the unspecified location case. 

Keywords: Circular statistics, Fisher information singularity, Rayleigh test for uniformity, 
skewed distributions, tests for symmetry 



1 Introduction 

Symmetry is a fundamental and ubiquitous structural assumption in statistics, underpinning 
most of the classical inferential methods, be it for univariate data on the real line or for circular 
data. Its acceptance generally simplifies the statistician's task, both in the elaboration of new 
theoretical tools and in the analysis of a given set of observations. For instance, the classical 
models for circular data, such as, e.g., the von Mises, cardioid, wrapped normal or wrapped 
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Cauchy distributions (see Mardia and Jupp 2000, Section 3.5) are all symmetric about their 
unique mode (this form of symmetry on the circle is called reflective symmetry). However, 
quoting Mardia (1972, p. 10), "symmetrical distributions on the circle are comparatively 
rare", and recent years have shown an increasing interest in non-symmetric models (see, e.g., 
Umbach and Jammalamadaka 2009, Kato and Jones 2010, Abe and Pewsey 2011 or Jones and 
Pewsey 2012); therefore, it is all the more important to be able to test whether the hypothesis 
of symmetry holds or not, in order to know whether the classical, well explored, or rather 
the modern, less explored, models shall be used. Since circular distributions are encountered 
in several domains of scientific investigation, with particular emphasis on the analysis of (i) 
phases of periodic phenomena (physics, biology, etc.) and (ii) directions (animal movements as 
a response to some stimulus, pigeon homing, earth sciences, etc.), practical examples needing 
tests for circular symmetry about a specific axis/direction associated with the experimental 
set-up under consideration are all but rare. 

On the real line, testing for symmetry about a fixed center (the median) is a classical yet 
timeless issue, and consequently there exists an abundance of such tests. Essentially, these 
tests can be distributed in two distinct categories. The first class contains nonparametric tests 
which are based on given characteristics of the null hypothesis of symmetry; famous examples 
are the Cramer- von Mises type tests of Rothman and Woodrofe (1972) or Hill and Rao (1977), 
the Kolmogorov-Smirnov type test of Butler (1969), the runs test of Mc Williams (1990) or its 
modified version by Modarres and Gastwirth (1996), to cite but these. While being able to detect 
asymmetry against (usually) any type of skew alternative, such tests suffer from nonparametric 
rates of convergence, hence need in general a large number of observations in order to become 
powerful. Tests belonging to the second class are not universally consistent but instead rather 
focus on some favored alternatives, against which they achieve (semi)parametric consistency 
rates and sometimes even are (semi)parametrically optimal; examples of this category are 
provided in Kozubowski and Panorska (2004), Cassart et al. (2008) or Ley and Paindaveine 
(2009). If one further considers shifts in location as "skew" alternatives, then the classical sign, 
Wilcoxon and van der Waerden tests (see Hajek and Sidak 1967), which respectively achieve 
optimality under double-exponential, logistic and normal distributions, also belong to this second 
category. 

In the circular case, the null hypothesis of symmetry is way less explored. There exist 
essentially three proposals for such tests: 

• Schach (1969) constructs locally optimal linear rank tests against rotation alternatives, 
the circular analogue of a linear shift alternative. His construction comprises the circular 
sign and Wilcoxon tests. 

• Universally consistent tests from the linear setting have been adapted to the circular case 
(such as the celebrated runs tests, see Pewsey 2004). 

• A "true" test for circular symmetry has been studied in Pewsey (2004) by having recourse 
to the second sine moment about the fixed median direction, a classical measure of circular 
skewness first proposed by Batschelet (1965). 

(One should not confuse the problem of testing for reflective symmetry treated here with that of 
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testing for /-fold symmetry on the circle; this issue has been addressed in Jupp and Spurr 1983 
(see also Mardia and Jupp 1999, page 146).) The scarcity of existing tests for circular symmetry 
might at first sight seem puzzling as one may be tempted to say that all tests for (linear) 
univariate symmetry should be adaptable to the circular setup (such as done for the runs tests) , 
replacing the real line by the compact [—tt, tt] . However, this translation from one setup to 
the other is not so straightforward, due to several facts including that the points at tt and —tt 
coincide. Furthermore, replacing the circle by the compact set [— tt, tt] requires the choice of an 
arbitrary "reference point" on the circle which will play the role of the zero on the real line; 
this sounds unrealistic. As stated in Pewsey (2004), when the observations are distributed on a 
large arc of the circle, it is likely that adapted tests suffer from a loss of power. It seems also 
very unlikely that optimal tests on the real line will keep their optimality features on the circle, 
as nothing a priori ensures that they behave well against the (certainly complicated) wrapped 
versions of the univariate skew distributions they were designed for. Thus, except for rotation 
alternatives, the category of optimal tests for reflective symmetry appears to be empty. 

In view of the emptiness of the second category of tests and the growing interest in skew 
circular distributions, our aim in the present paper is to fill in this gap by proposing tests for 
circular reflective symmetry about a fixed center that behave extremely well against a certain 
(general) type of skew alternatives. More precisely, we shall build locally and asymptotically 
optimal (in the maximin sense) tests for symmetry against k -sine- skewed alternatives (Umbach 
and Jammalamadaka 2009, Abe and Pewsey 2011), a broad class of recently proposed skew 
circular distributions that has received an increasing interest over the past few years (see 
Section 2 for a description). In a nutshell, these skew distributions are obtained by perturbation 
of a base symmetric distribution via a factor involving sines and a parameter to regulate 
skewness. Apart from the general interest in these skew circular distributions, the motivations 
for this choice are mainly twofold. First, they are the circular analogues of the most studied 
and most used skew distributions on the real line, namely the skew-symmetric distributions (see 
Azzalini and Capitano 2003 or Wang et al. 2004) inspired from the skew-normal distribution 
proposed in the seminal paper Azzalini (1985). Second, as we shall see, the resulting test 
statistics are based on the trigonometric sine moments, hence provide this classical measure of 
circular skewness as well as the test of Pewsey (2004) with so far not known optimality properties. 
As nice by-products, our findings also enable us to (i) build optimal tests for uniformity against 
cardioid alternatives, and (ii) discuss Fisher singularity issues exactly as in the linear case. 

The backbone of our approach is the Le Cam methodology which, although of linear nature, 
lends itself well for a transcription to circular settings (and even, with much more complications, 
to data living on unit hyperspheres in higher dimensions, see Ley et al. 2013). In a first stage, we 
will obtain optimal parametric tests, and then, by means of studentization arguments, we shall 
turn them into semiparametric ones, valid under the entire null hypothesis of symmetry and 
optimal not only, as is usually the case, under the symmetric base distribution their parametric 
antecedents are based on, but uniformly optimal under any given symmetric base distribution. 
We will hence derive, as Schach (1969), a family of fully efficient semiparametric tests which, 
in our case, are always optimal. For a given density, our tests will thus behave asymptotically 
as the likelihood ratio tests, but they clearly improve on the latter by their simplicity and the 
fact that, thanks to the Le Cam approach, one can derive explicit power expressions against 
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sequences of contiguous skew alternatives. 

The paper is organized as follows. In Section 2, we first describe the family of k -sine-skewed 
distributions, then establish their ULAN property in the vicinity of symmetry, the crucial step 
in the Le Cam approach, and discuss some aspects of this property. In Section 3, we construct 
our optimal tests for reflective symmetry about a known center and investigate their asymptotic 
properties. The finite-sample performances of all our tests for reflective symmetry are evaluated 
and compared to existing tests in a large Monte Carlo simulation study, see Section 4. Our 
tests are then applied on a famous real-data set in Section 5. In Section 6, we show how our 
findings also allow us to produce tests of uniformity against cardioid alternatives. The Fisher 
information singularity issue is tackled in Section 7. We conclude the paper with final comments 
and an outlook on the case of non-specified median direction in Section 8, and an Appendix 
collects the technical proofs. 

2 k -sine-skewed distributions and the ULAN property 

2.1 k -sine-skewed densities 

As briefly depicted in the Introduction, A; -sine-skewed distributions are obtained by perturbation 
of a base symmetric density. Define the collection 

T := Uo:fo(x)>0eL.e.,f Q (x + 2irk) = f Q (x)VkeZ,f (-x) = fo(x), 
/o unimodal at 0, J fo(x)dx = ll 

of unimodal reflectively symmetric (about the zero direction) circular densities. The most 
well-known representatives of this collection are the von Mises, cardioid or wrapped Cauchy 
distributions, with respective densities fyM K (x) '■= 2-kI (k) ex P( /t cos(x)) for k > (Jo stands 
for the modified Bessel function of the first kind and order zero), fcA e (x) '■= ^(1 + ^cos(x)) for 
t E (0,1), and fwc p (x) := 1+p i_2 pcos ( x ) for P G (0, 1). A location parameter 9 G [-7r,7r] is 
readily introduced as center of symmetry, leading to densities f(x — 6),x G [— vr,7r], with mode 
9 . Inspired by the classical one-dimensional skewing method of Azzalini and Capitanio (2003), 
Umbach and Jammalamadaka (2009) have skewed such symmetric densities /o by turning them 
into 

2f (x-9)G(oj(x-9)), xe[-7r,7r], 

where G(x) = J_ g(y)dy is the cumulative distribution function (cdf) of some circular 
symmetric density g and u is a weighting function satisfying for all x G [— ir, ir] the three 
conditions lo{— x) = —uj(x), uj(x + 2irk) = u;(x)yk G Z, and |w(x)| < ir . This construction being 
too general and for the sake of mathematical tractability, Umbach and Jammalamadaka have 
particularized their choice to G(x) = (tt + x)/(2ir) , the cdf of the uniform circular distribution, 
and uo{x) = A7rsin(fca;), k G No, with A G (—1, 1) playing the role of a skewness parameter. This 
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finally yields what we call the k -sine-skewed densities 

x^f^ x (x):=f (x-6)(l + \sm(k(x-e))), x€[-7r,7r], (2.1) 

with location parameter 9 G [— tt,it] and skewness parameter A G (—1,1). When A = 0, no 
perturbation occurs and we retrieve the base symmetric density, otherwise (2.1) is skewed to the 
left (A > 0) or to the right (A < 0). Further appealing properties of k -sine-skewed distributions 
are that fg X (6 — x) = fg_ x (6 + x) , fjf x (9) = /o(0) independently of the value of A, and the 
two endpoints, fg X (9 — n) and fg X {9 + vr), coincide. However, for k > 2, fg X happens to be 
multimodal, whereas, for k = 1, multimodality only rarely occurs. This explains why Abe and 
Pewsey (2011) have restricted their attention to the study of the densities 

x^f e ,x(x):=fo(x-e)(l + \sm(x-6)), x€[-7r,7r], (2.2) 

which they have called sine-skewed circular densities (hence our terminology A; -sine-skewed 
densities for general k). Abe and Pewsey have shown in their paper under which conditions the 
densities (2.2) happen to be multimodal. In the present paper, we establish all our theoretical 
results and propose tests for general k -sine-skewed distributions. Note that, when /o is the 
circular uniform density, then (2.2) corresponds to the cardioid density /ca a with mode at 
9 + 7r/2(mod 2ir) , hence, in passing, we will as well construct and analyze an optimal test for 
uniformity against the cardioid distribution. 

Sine-skewed (and A: -sine-skewed) distributions lend themselves pretty well for modeling real- 
data phenomena. Aside from Abe and Pewsey (2011) where this aspect is thoroughly described, 
these skew-circular distributions have been used, inter alia, in the analysis of the CO2 daily 
cycle in the low atmosphere at a rural site (Perez et al. 2012) and of forest disturbance regimes 
(Abe et al. 2012). This, combined with the motivations stated in the Introduction, makes k- 
sine-skewed distributions a natural choice as asymmetric alternatives in the construction of tests 
for circular reflective symmetry. 

2.2 The ULAN property for A; -sine-skewed densities 

As explained in the Introduction, we shall have recourse to the Le Cam methodology in order 
to construct locally and asymptotically optimal tests for reflective symmetry against k -sine- 
skewed alternatives. For the sake of generality and in view of future research (see Section 8), 
we here do not assume 9 to be fixed. This of course contains the #-fixed case, which we need 
in this paper. Let X±, . . . ,X n be i.i.d. circular observations with common density (2.1) (that 
is, the Xi's are angles). For any symmetric base density /o G T and any k G No, denote by 
F'tf- jb k , where 1? := (9, A)' G [— 7r, tt] x (—1, 1) , the joint distribution of the n-tuple X\, . . . , X n . 
Since, for A = 0, the density fg X reduces to /o and hence does not depend on k, we drop the 

index k and simply write P^™| at 1? = t?o := (9,0)' . Any pair (fo,k) induces the parametric 
location-skewness model 

*ft:= h^-] *(-M)}, 
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whereas any k G No induces the nonparametric location-skewness model := Uf oe jrVj^\. 

The very first step of our construction of tests for symmetry about a fixed center consists 
in establishing the Uniform Local Asymptotic Normality (ULAN) property, in the vicinity of 
symmetry (i.e., at A = 0), of the parametric model rf^ k - This property of A; -sine-skewed 
distributions happens to be interesting per se, as it paves the way to numerous other applications 
of the Le Cam theory (such as, e.g., the construction of tests for symmetry about an unspecified 
center or of the celebrated one-step optimal estimators, see e.g. van der Vaart 2002). ULAN 
requires the following mild regularity condition on the base densities fo . 

Assumption (A). The mapping x i— > fo(x) is a.e.-C 1 over [— tt, tt] (or equivalently over R by 
periodicity) with a.e.-derivative fo . 

Most classical reflectively symmetric densities satisfy this requirement. Note that the C 1 
condition over a compact combined with the fact that fo > entails that, letting ipj = — fo/ fo, 
the Fisher information quantity for location If := f^ip 2 o (x)fo(x)dx is finite. ULAN of the 

parametric model ~P^\ with respect to "& = (9, A)', in the vicinity of symmetry, then takes the 
following form. 

Theorem 2.1. Let fo £ J 7 and k G No, and assume that Assumption (A) holds. Then, for 
any 9 G [— tt,tt], the parametric family of densities Vj?\ is ULAN at i?o = {^,0)' with central 
sequence 



A (n) {9) ._ | foM 
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V™ \ s'm(k(Xi 
and corresponding Fisher information matrix 

T fo,k ■ = 

where r /0)fe;11 := I fo , ^f ,k;l2 ■= ~ f- w sin(kx)f (x)dx and 

r/ ,fc;22 := fl~ sm 2 (kx)fo(x)dx . More precisely, for any = 9 + 0{n~ l l 2 ) and for 

any bounded sequence = (r| n , ri n )' G M 2 such that n~^l 2 T^ belongs to (—1,1) 

and #( n ) + n~ l / 2 T^ remains in [— tt,tt] (with a slight abuse of notation since, e.g., 
—tt — e = tt — 6 for e > 0; what we mean is that the perturbation of the circu- 
lar location parameter 9 is such that it still lies on the unit circle), we have, letting 

A (n) = T (»)' A W fc (0(»O) _ (i/2) r W'r /o , fc rW + op(l) 

(2.3) 
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and AfJ k (9( n ))AX 2 (p,T f0tk ), both under P$ n))0) , /0)& as oo. 

The proof is given in the Appendix. One easily sees that the Fisher information for skewness 
r/ ,fc;22) and hence the cross-information quantity r / 0i fc ; i2, is finite by bounding sin 2 by 1 under 
the integral sign. Note that the constant k has no effect on the validity of Theorem 2.1. Note also 
that A^"^) does not depend on / , a fact that will become of great interest in the sequel. With 
this ULAN property in hand, we are ready to derive our optimal tests for reflective symmetry 
about a fixed center 9, as explained below in Section 2.3. Moreover, since we do not fix 9 in 
Theorem 2.1, our result also paves the way for deriving optimal tests for symmetry about an 
unknown center; see Section 8. 

We conclude the present section with a brief discussion on the minimal conditions required 
for having the ULAN property. Indeed, in view of the proof of Lemma A.l in the Appendix 
which is the main step to demonstrate Theorem 2.1, Assumption (A) can be further weakened 
into 

1/2 

Assumption (A mm ). The mapping x i— >■ f (x) is differentiable in quadratic mean over 

1 /2 

[— 7r, 7r] (or equivalently over R by periodicity) with quadratic mean or weak derivative (/ )'(x) 

1/2 1/2 

and, letting vpf (x) = — 2(/ )'{%)/ fo ( x ), the Fisher information quantity for location 
Jfo '■= fl n ^j { x )fo{x)dx is finite. 

1/2 

Quadratic mean differentiability of / , a classical requirement in the Le Cam framework, 
means that _(/q (x — h) — f^ 2 (x) — htpf (x)) 2 dx = o(h 2 ) as h — > 0, which corresponds 

1 /2 

exactly to the integral in expression (A. 9) of the proof of Theorem 2.1 with (/ )' instead 

1 1/2 

of 2/0(^)7/0 i x ) an d hence is the minimal condition in order to have the ULAN property of 
the parametric model Vj?\- Note that, under Assumption (A), these two derivatives of course 
coincide well as tfjf and (pf , and Jf = If (as already mentioned above, the C 

condition ensures finiteness of If , while in the weaker Assumption (A mm ) one needs to ask 
that Jf < 00). 

2.3 Constructing Le Cam optimal tests on basis of the ULAN property 

The central idea of the Le Cam theory we are using here is the concept of convergence of 
statistical models (experiments in the Le Cam vocabulary). Quoting Le Cam (1960), "the 
family of probability measures under study can be approximated very closely by a family of 
a simpler nature". The ULAN property is an essential ingredient in this approximation, as it 
allows to deduce that (see Le Cam 1986 for details) our parametric location-skewness model 
is locally (around (#,0)') and asymptotically (for large sample sizes) equivalent to a 
simple Gaussian shift model. Intuitively, this follows from the fact that the likelihood ratio 
expansion (2.3), if we neglect the remainder terms, strongly resembles the likelihood ratio of a 
Gaussian shift model A/2(r j ^T^ n \T j ^) with a single observation denoted by A^ fc . Since the 
optimal procedures for Gaussian shift experiments are well-known, we can translate them into 
our circular location-skewness model and hence obtain inferential procedures that happen to be 
asymptotically optimal (here, in the maximin sense). 

In the present paper, we shall employ this working scheme for testing the null hypothesis 
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T~Lq of symmetry about a known central direction 6 G [— 7r, 7r]. As explained above, our 
procedures will be (asymptotically) optimal against a fixed k -sine-skewed alternative (2.1). 
For the considered testing problem, we first construct fo -parametric tests for %o-/o = $a)'-fo : 
the optimality of these tests under the base density /o is thwarted by the fact that they are 
only valid (in the sense that they meet the asymptotic level- a constraint) under /q. In order 
to palliate this non- validity outside /o , we have recourse to a classical studentization argument 
allowing us to turn our parametric tests into tests for the semi- parametric null hypothesis 
%q = U/- oe jrP|"^ y . The next section contains the detailed derivations of these tests. 

3 The test statistic and its asymptotic properties 

Fix 9 G [— vr,7r]. The /o-parametric test <f>^.f for circular reflective symmetry about a known 
central direction 6 we propose rejects %o-/o a * as y m Pt°ti c level a whenever the statistic 

n (n);e _ l A SWI _ |n- 1/2 Eg=l«m(fe(^-g))| d) 
^fo;k — 1/2 ~ p l/2 V A > 

exceeds z a / 2 , the a /2 upper quantile of the standard normal distribution (tests for reflective 
symmetry against one-sided alternatives of the form A > or A < are built similarly). It 
follows from the Le Cam theory that this test is locally and asymptotically maximin for testing 
the null %o-,fo a § arns t ^i-f k := ^A^0e(-i,i)P(e a)'-/ k (this optimality does not hold for k' ^ k). 

Now consider go G F. Under P(0o)'-g o ' k-2$) ls asymptotically normal with mean and 
variance T ^fo^k-pi- It is therefore natural to consider the studentized test <p k that 

rejects (at asymptotic level a) the null of circular reflective symmetry as soon as 

Q *(ny,e . = \n- 1 / 2 E7=iMHx i -e))\ 

k ' (n~i£r=iSm 2 (A^-*))) 1/2 

exceeds z a /2 ■ We attract the reader's attention on the fact that this very simple test statistic 
does no more depend on /n (hence the omitment of the index /n in <^t ). This is of course 
due to the fact that the central sequence for skewness, a!"]^), does not depend on /q. This 
remarkable fact implies that all parametric tests $J]y£,k £ N, lead to the same studentized 

test statistic (j) k , which therefore inherits optimality from its parametric antecedents under 
any base symmetric distribution! This nice property as well as the asymptotics of such tests, 
under any /q G T , follow from the ULAN property in Theorem 2.1 and are summarized in the 
following result (see the Appendix for a proof). 

Theorem 3.1. Let k G N . Then, 

(i) under U/ gjrP^L, , Q* k n ^' 9 -5- Af(0, 1) as n — >■ oo, so i/iai £/ie tesi (f)* k n ^' e has asymptotic 
level a under the same hypothesis; 

fii) under , ,„•, with fn G J- and k! G No, Qt^' 9 is asymptotically normal 

iuif/i mean Tj^^C/o (/c, A; / )t2 and variance 1, where t 2 = liuin-yoo t 2 and Cf (k,k') : = 



8 




Figure 1: Power curves, as a function of t 2 , of the studentized test 4> k for k = 2 against 

(n) 

local alternatives P , for fn the von Mises density with concentration parameter 

1 and for k' equal to 1 (blue line), 2 (red line) and 3 (yellow line). 

sin(kx) sin(k' x) fo(x) dx (which is finite); 
(Hi) for all /o 6 T , Q*^ ,e = Qf^.f + °p(l) as n — >■ 00 under P(g o)' /o ' so ^ ia ^ ^ e s ^ u ^ en ^' iZe - ( ^ 
test ^)*S n '' S is locally and asymptotically maximin, at asymptotic level a, when testing H® against 
alternatives of the form U A ^ e(-i,i) u / 6^ P (^A)';/ ,fe ■ 

Theorem 3.1 (i) shows that indeed the studentized test 0t is valid under the entire null 
hypothesis H®, hence is asymptotically distribution-free. Note the uniform (in /o, not in k) 
optimality of our studentized test. 

For the sake of generality, we have also considered above alternatives where the natural k is 
replaced by some k! possibly different from the k used in the construction of our tests. Point (ii) 

allows us to give the explicit power of &*S n ^' e against the local alternatives P^ , , , : 

k (S,n- 1 /Vf ) )' i /o,it' 

1 - $ (*«/2 " (r/ ,fc;22)- 1/2 C /0 (fc, k')r 2 ) + d> {-z a/2 - (T fo>k]2 2)- 1/2 C fo (k, k')r 2 ) , 

where $ stands for the cdf of the standard Gaussian distribution. For the sake of illustration 
we have plotted, in Figure 1, the power as a function of t 2 for /o the von Mises density with 
concentration parameter 1 and for k = 2 and k' = 1,2,3. Quite interestingly, when k = 2, the 
statistic cf)*^ 1,9 coincides exactly with the so-called "b2-star" test proposed in Pewsey (2004). 
We have thus shown that this test enjoys maximin optimality features against 2 -sine-skewed 
alternatives, and provided its asymptotic powers against contiguous alternatives. This not only 
complements, but also gives further insight into the b2-star test. Finally, the very simple tests 
we have obtained are also easy to interpret as they are based on trigonometric sine moments, 
which are classical measures of skewness for circular data (see, e.g., Batschelet 1965). 
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4 Monte Carlo simulation study 



The aim of this section consists in examining the finite-sample performances of the proposed 
testing procedures for reflective symmetry (check of the nominal level constraint under distinct 
forms of reflective symmetry and investigation of power properties under various forms of 
asymmetry). To this end, we have generated N = 3,000 independent samples of size n = 100 
from reflectively symmetric and increasingly skewed (A > 0) circular distributions, and run our 
tests as well as a natural competitor of the literature (the modified runs test of Pewsey 2004) 
under two-sided form at the asymptotic level a = 5% . 

Let us fix the center of symmetry 9 to 0. We have performed our tests 0*^'°, i^"^' 
and cf)^^' as well as the modified runs test 0^ drun with p = 0.6. The lat- 
ter test, initially proposed on the real line by Modarres and Gastwirth (1996) as 
a modification of the classical runs test of McWilliams (1990), has been adapted in 
Pewsey (2004) to the circular case (our choice of p = 0.6 also stems from that pa- 
per). The reason for choosing (^odrun as competitor lies in the simulation study driven 
in Pewsey (2004). Indeed, recall that cf)^^' coincides with the so-called b2star test of 
Pewsey (2004); in that paper, the author has compared the performances of the test </>2^' 

(n) 

against various competitors, and the conclusion was that only <^rnodrun could compete with 
02^'°, whence our restriction. Of course, we here choose other settings as in Pewsey (2004), 
and hence also add new simulation-based information on Pewsey's b2star test. We also remark 
that we consider here k = 1, 2, 3 for our tests because these values are able to capture both skew 
unimodality (k = 1) and multimodality, but do not lead to too many oscillations of the sines 
within [— 7r, 7r] . 

As reflectively symmetric distributions representing the null hypothesis, we have considered 
the von Mises laws /vMi and /vm 10 , the cardioid /ca 5 > the wrapped Cauchy /wc 5 as wen 
as a mixture of two /vm x and two /vm 10 von Mises laws with, in each case, respective centers 
at — 7r/4 and 7r/4 and mixing probability 0.5. The latter mixture is used in order to assess 
the performances of our tests under bimodality. The densities /vMi and /vMi have then been 
turned into their l-,2- and 3-sine-skewed versions, whereas fcA 05 an d /wc 05 have become 

1- and 2 -sine-skewed. More concretely, in each case the skewness parameter A increases from 
zero to successively positive values. The bimodal mixture of von Mises laws has been skewed by 
simply shifting the center 7r/4 to 7r/4+A. In order to also investigate other forms of perturbation 
of symmetry, we have applied the Moebius transform of Kato and Jones (2012) to /vMi an d 
/vMio with r = 0.5 (in short, it turns X into A + 2arctan (oj r tan ( ^f^ )) with u r = jz^p)- 
The rejection frequencies are reported in Table 1 for 1 -sine-skewed alternatives, Table 2 for 

2 - sine-skewed alternatives and Table 3 for Moebius, von Mises mixtures and 3-sine-skewed 
alternatives. 

All four tests meet the 5% nominal level constraint under each reflectively symmetric density 
considered, even under bimodality, and seem to be unbiased. Under k -sine-skewed alternatives, 
the theoretical optimality features of our tests <f> k are confirmed. Quite remarkably, when 
the observations are highly concentrated (/vMi case), the differences between the three tests 
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Table 1: Rejection frequencies (out of N = 3,000 replications), under various reflectively 
symmetric and 1 -sine-skewed distributions, of the optimal tests (f>*^'° . '° and (j)^ n ^'° as 
well as of the modified runs test with p = 0.6. The tests are performed at level a = 5%. 



1-sine-skewed /vmi 



Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0517 


0.2633 


0.7600 


0.9873 




0.0567 


0.0927 


0.2277 


0.4373 


,*(n);0 


0.0497 


0.0537 


0.0593 


0.0713 


ai (n) 

t modrun 


0.0577 


0.0747 


0.1620 


0.3667 


1-sine-skewed /vm 10 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 


,.(n);0 
<Pl 


0.0460 


0.0993 


0.2333 


0.4573 


^(„) i0 


0.0470 


0.0963 


0.2303 


0.4523 




0.0507 


0.0900 


0.2137 


0.4203 


<A (n) 

t modrun 


0.0583 


0.0610 


0.0880 


0.1387 


1-sine-skewed fcA 5 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0493 


0.2993 


0.8220 


0.9960 




0.0497 


0.0620 


0.0940 


0.1883 




0.0497 


0.0447 


0.0527 


0.0460 


(A (n) 

t modrun 


0.0487 


0.0630 


0.1243 


0.2680 


1-sine-skewed /wco.s 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0510 


0.2417 


0.6943 


0.9703 


^(„) i0 


0.0507 


0.0853 


0.1973 


0.3827 




0.0513 


0.0637 


0.0860 


0.1250 


r modiu n 


0.0523 


0.0683 


0.1493 


0.3410 
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Table 2: Rejection frequencies (out of N = 3,000 replications), under various reflectively 
symmetric and 2 -sine-skewed distributions, of the optimal tests (f>*^'° . '° and (j)^ n ^'° as 
well as of the modified runs test with p = 0.6. The tests are performed at level a = 5%. 



2-sine-skewed /vm 1 



Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0460 


0.1100 


0.2423 


0.4917 




0.0473 


0.3040 


0.8287 


0.9950 


,*(n);0 


0.0550 


0.0907 


0.2427 


0.4823 


ai (n) 

t modrun 


0.0577 


0.0677 


0.1137 


0.2487 


2-sine-skewed /vm 10 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 


,.(n);0 
<Pl 


0.0450 


0.1813 


0.5753 


0.9133 


^(„) i0 


0.0457 


0.1833 


0.5843 


0.9157 




0.0467 


0.1800 


0.5653 


0.9013 


t modrun 


0.0533 


0.0627 


0.1730 


0.4130 


2-sine-skewed fcA 5 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 


£(„),<> 


0.0507 


0.0633 


0.1027 


0.1840 




0.0507 


0.2717 


0.8313 


0.9917 




0.0520 


0.0700 


0.1097 


0.1787 


(A (n) 

t modrun 


0.0553 


0.0703 


0.1350 


0.2853 


2-sine-skewed /wc .5 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0493 


0.0920 


0.2347 


0.4543 


^(„) i0 


0.0597 


0.2710 


0.7913 


0.9913 




0.0563 


0.0937 


0.2600 


0.5230 


r modiu n 


0.0480 


0.0640 


0.1367 


0.2587 



12 



Table 3: Rejection frequencies (out of N = 3,000 replications), under various reflectively 
symmetric and various skewed distributions, of the optimal tests (j>*^'°, (p^ n ^'° and (p^ n ^'° as 
well as of the modified runs test ^Idrun w hh p = 0.6. The tests are performed at level a = 5%. 



Moebius transformed /vmi 



Test 


A = 


A = 0.2/3 


A = 0.4/3 


A = 0.2 


,*{n);0 
<Pl 


0.0447 


0.1263 


0.3437 


0.6487 


^(n);0 


0.0480 


0.1567 


0.4603 


0.7923 


^(n);0 


0.0497 


0.1637 


0.4567 


0.7720 


"modrun 


0.0563 


0.0570 


0.0777 


0.1133 






Moebius transformed /vm 


10 


Test 


A = 


A = 0.02 


A = 0.04 


A = 0.06 


^(n);0 


0.0510 


0.2303 


0.6553 


0.9357 


^(n),0 


0.0517 


0.2327 


0.6563 


0.9363 


^) i0 


0.0520 


0.2337 


0.6613 


0.9380 




0.0527 


0.0850 


0.2020 


0.4623 


Skewed /m, mixtures 


Test 


A = 


A = 0.4 


A = 0.8 


A = 1.2 


<f>l 


0.0470 


0.1007 


0.1460 


0.0980 




0.0457 


0.0550 


0.1240 


0.2550 




0.0493 


0.0470 


0.0537 


0.0543 


<A Cn) 

r modrun 


0.0453 


0.0600 


0.0787 


0.0890 






Skewed /vm 


10 mixtures 




Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0507 


0.1303 


0.2463 


0.3560 




0.0510 


0.0693 


0.3500 


0.9573 


^(n);0 


0.0490 


0.6607 


1.0000 


1.0000 


rmodrun 


0.0553 


0.2373 


0.8027 


0.9923 


3-sine-skewed /vmi 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0563 


0.0450 


0.0603 


0.0787 




0.0483 


0.0963 


0.2420 


0.4657 


^(»);0 


0.0503 


0.3027 


0.8247 


0.9920 


< modrun 


0.0550 


0.0657 


0.1343 


0.2840 


3-sine-skewed /vm 10 


Test 


A = 


A = 0.2 


A = 0.4 


A = 0.6 




0.0480 


0.2297 


0.6963 


0.9687 




0.0470 


0.2410 


0.7347 


0.9783 




0.0473 


0.2520 


0.7530 


0.9827 


<t> {n \ 

" modrun 


0.0557 


0.0890 


0.1947 


0.5207 
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vanish, whereas under certain k' -sine-skewed densities the test 4>* k for k ^ k! exhibits low 
powers (especially when combining the indices 1 and 3). We moreover notice that all our tests 
are powerful under the Moebius transformed skew densities and even under skewed von Mises 
mixture distributions with high concentration parameter k , which shows that the proposed tests 
not only act well under the laws they are designed for. As an overall summary, we see that our 
three tests generally outperform the modified runs test. 

5 A real-data application 

In this section, we apply our optimal tests for reflective symmetry on a well-known data set in 
animal orientation problems. This data set stems from an experiment with 730 red wood ants 
(Formica rufa L.) described in Jander (1957). Each ant was individually placed in the center 
of an arena with a black target positioned at an angle of 180° from the zero direction, and the 
initial direction in which each ant moved upon release was recorded to the nearest 10° . Thus it 
is clear that the experimental design suggests a natural median direction, a fact that is clearly 
corroborated by the graphical representation of the data in Figure 2. The natural question of 
interest is whether the directions chosen by the ants are symmerically distributed around the 
median direction representing the black target. This data set has been discussed, inter alia, in 
Pewsey (2004), Umbach and Jammalamadaka (2009) and Abe and Pewsey (2011). 

By the experimental design, this real-data set happens to be a very good candidate for testing 
circular symmetry about a known median direction (assuming the median direction unknown 
might even be inappropriate here and the corresponding tests for symmetry will not be as 
powerful as tests for symmetry about a fixed direction). In view of the data plot in Figure 2, we 
see that the underlying density might be multimodal rather than unimodal, indicating that the 
tests (j)^ 12 ^' and (j)^ n '' might be more powerful in the present situation than </>*^'° (we refer to 
Abe and Pewsey 2011 for a discussion on the conditions under which 1 -sine-skewed distributions 
are unimodal or multimodal). Indeed, c/>*^'° yields a p-value of 0.7781, while (j)^^' and 0g^'° 
respectively give p-values 0.0107 and 0.0131. This clearly provides evidence that the data are 
in fact not symmetrically distributed around the median direction of 180° . Pewsey (2004) has 
obtained the same conclusion with his b2star test, but our conjunction of the tests </> 2 and 
^*(n),o p rov ^ eg f ur ther information and evidence. In Abe and Pewsey (2011), the authors 
notice that neither the symmetric nor the 1-sine-skewed distributions they have considered 
provide an adequate fit to this data. Their findings are not a surprise: according to (p-^ ' , 
1-sine-skewed densities are not preferable over symmetric ones, while our other tests reject the 
hypothesis of reflective symmetry up to the 2% level. This shows that, most probably, the ant 
data are best fitted by 2- or 3-sine-skewed distributions. 

6 Optimal test for uniformity against cardioid alternatives 

In this section, we show that the parametric test statistic (3.4) not only serves for efficient 
testing of reflective symmetry but also for testing the null hypothesis of uniformity against 

cardioid alternatives about a fixed central direction. Mardia and Jupp (2000) explain that there 
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Figure 2: Raw circular plot of the Jander (1957) data set recorded during an orientation 
experiment with 730 red wood ants. Each dot represents the direction chosen by five ants. 

exist circumstances under which one wishes to test uniformity against a certain fixed direction; 
such situations are the motivation behind our optimal test for uniformity. As described in 
Section 2.1, for fo(x) = the uniform density, (2.2) with k = 1 corresponds to the cardioid 
density with mode at 9 + tt/2 € [— 7r,7r]. Of course, we cannot directly use our ULAN result 
for this hypothesis, as the uniform distribution does not favour one particular direction, hence 
9 would not be identified under %Q mf (the uniform density does not belong to J-). However, 
writing Pj^jjf for the joint distribution of an n-tuple of uniform observations X±, . . . ,X n and 
(o \)' ^ or ^ ne j° m t distribution under the alternative T-L^ Ax ' 9 of cardioid distribution, and using 
the Taylor expansion log(l + x) = x — ^ + o(x 2 ) around 0, we obtain that 

n 

log(l + ra~ 1/2 r 2 sm{Xi - 9)) 

i=l 

n 2 

r 2™~ 1/2 S sin ( Xi " e ^ ~ y E [ sm2 (^i - *)] + Mi) 
i=i 

as n — > oo under Pj^jL, which is nothing but a LAN decomposition. Hence, for fixed 9, an 
optimal test for testing 7-Lq : A = against T~L 1 A ' : A 7^ shall be based on 

n 

1=1 

that is, on (3.4) with k = 1 and r / unifi fc ; 22 = \ - The results of Theorem 3.1 remain valid under 
this special parametric setup here (with, of course, the required modifications such as restricting 
t° fo = /unif m Parts (i) and (iii) and restricting k' to 1 and /o to the cardioid density in 
Part (ii)), and therefore provide the asymptotic behavior of Q^tL.g and the working mechanism 

of the associated test of uniformity 0™jf.g (which is a one-sided test in view of the definition 



( dF (e,n-y 2 T 2 y\ 
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of the cardioid density). Now, recall that the central direction is not 9 but 9 + vr/2, hence the 
statistic (6.6) should rather read 

n n 

Qulw = yftn-W^saiiXi -(9 + vr/2) + vr/2) = v^rT^^cospQ - (9 + vr/2)), 

i=l i=l 

which corresponds exactly to the expression of the Rayleigh (1919) test of uniformity against 
a fixed central direction 9 + 7r/2, which is originally constructed to be efficient against von 
Mises alternatives, see Mardia and Jupp (2000), pages 98 and 99. Thus, our findings now 
further indicate that the Rayleigh test is also optimal against cardioid alternatives. Note that 
Theorem 3.1 provides its asymptotic null and non-null behavior. 



7 Singularity of the location-skewness Fisher information ma- 
trix 

Besides its numerous favourable properties, the skew-normal distribution of Azzalini (1985) is 
also sadly famous for having a singular Fisher information matrix in the vicinity of symmetry, 
due to the collinearity of the scores for location and skewness. A vast literature has been 
devoted to the analysis of the reasons for this singularity, to possible cures and to the study of 
which other skew-symmetric distributions suffer from the same drawback. Indeed, a singular 
information matrix violates the assumptions for standard Gaussian asymptotics of the maximum 
likelihood estimators and precludes, at first sight, any nontrivial test of the null hypothesis of 
symmetry. Among the papers having dealt with this singularity issue are Azzalini (1985), 
Azzalini and Capitanio (1999), Pewsey (2000), DiCiccio and Monti (2004), Chiogna (2005), 
Azzalini and Genton (2008), Ley and Paindaveine (2010) and Hallin and Ley (2012). The latter 
paper exactly determines those skew-symmetric distributions that suffer from Fisher singularity. 

The present section can be inscribed into this stream of literature, as it discusses and solves 
the same problem for A; -sine-skewed circular distributions. Moreover, our investigation shall 
become very important when one wants to construct optimal tests about an unknown center 9 , 
as will be briefly discussed in the final section. Now, recall that the information matrix in the 
vicinity of symmetry is given by 

r = ( n n ^%( x )fo( x ) dx n n sm(kx)ip fo (x)f (x)dx \ 
h,k \ fl n sm(kx)ip fo (x)f (x)dx f* n sin 2 (kx)f (x)dx J ' 

This matrix is singular if and only if 

l P%( x )fo(x)dx\ (^j sin 2 (kx)f (x)dx^j = (J sm(kx)<pf (x)f Q (x)dx\ . (7.7) 

The Cauchy-Schwarz inequality readily yields that the equality sign " = " in (7.7) can be replaced 
by ">" with equality holding if and only if (ff (x) = asin(kx) for some real constant a. The 
latter easy-to-solve first-order differential equation then shows that an information singularity 
can only occur for base symmetric densities /o of the form cexp(| cos(fex)) for a £ R and 
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c > a normalizing constant. Now, bare in mind that the class of base densities J- we consider 
contains the condition of unimodality on /q , which directly rules out all values k > 2 and forces 
a to be positive. Hence, the only base symmetric density for which the Fisher information 
matrix Tf 0i ji is singular corresponds to fo(x) = cexp(Kcos(a;)) with k = % > a concentration 
parameter, hence to the famous von Mises circular density. We formalize this result in the 
following proposition. 

Proposition 7.1. Let /q be a symmetric base density belonging to T and satisfying Assump- 
tion (A), and consider k -sine- skewed densities of the form fo(x—9)(l+\sm(k(x—9))) . Then the 
Fisher information matrix associated with the parameters 9 € [— ir, ir] and A £ (-1,1) is singular 
in the vicinity of symmetry (that is, at A = ) if and only if k = 1 and fo(x) = cexp(«;cos(x)) 
with k > a concentration parameter and c > the normalizing constant, that is, if and only 
if one is considering sine-skewed von Mises densities. 

Taking a closer look, this singularity only under the prominent von Mises densities is not 
really surprising. Indeed recall that the construction of asymmetric circular distributions by 
Umbach and Jammalamadaka (2009) is inspired from the classical linear setup initiated by 
Azzalini (1985). For the univariate setting, it has been shown in Hallin and Ley (2012) that the 
information matrix of skew-symmetric densities of the form 2f{x)G{\x) (with self-explaining 
notations) is singular if and only if the base symmetric density / is Gaussian. The circular 
equivalent of the linear first moment x being sin(x) (since cos(:r) is not odd), it is intuitively 
reasonable that the circular equivalent of the linear Gaussian density, that is, the von Mises 
distribution, suffers from Fisher singularity in the vicinity of symmetry, which allows to better 
understand the Fisher singularity result encountered here. 

8 Final comments 

In this paper we have tackled the problem of testing circular reflective symmetry about a 
specified center. The tests we propose are uniformly (over the null hypothesis) locally and 
asymptotically maximin against A: -sine-skewed alternatives, asymptotically distribution-free and 
moreover of a very simple form. They furthermore exhibit nice finite sample behaviors. Now, 
as already mentioned before, it would also be of interest to adapt our procedures to the case 
of an unspecified center, and our general ULAN property provides the required theoretical 
background for constructing such tests. The crucial difference to the tests of the present paper, 
of course, lies in the fact that we will need to replace the unknown location 9 with an estimator 
9. Would the information matrix Tf ± be diagonal, then the substitution of 9 for 9 would 
have no influence, asymptotically, on the behavior of the central sequence for skewness (9) . 
However, the covariance T f k;i2 only rarely equals zero, hence a local perturbation of 9 has the 
same asymptotic impact on (9) as a local perturbation of A = . It follows that the cost of 
not knowing the actual value of the location 9 is strictly positive when performing inference on 
A; the stronger the correlation between 9 and A, the larger that cost. The worst case occurs 
of course when the information matrix is singular (see Section 7), which leads to asymptotic 
local powers equal to the nominal level a; more precisely, this situation entails that the best 
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possible test is the trivial test, that is, the test discarding the observations and rejecting the 
null of reflective symmetry at level a whenever an auxiliary Bernoulli variable with parameter 
a takes value one. 

Now, in order to take into account the aforementioned cost of not knowing 9 , one can replace 
the central sequence A^(9) with the so-called efficient central sequence 

\( n ) e fffa\ A( n Vm r/o,fc;i2 A (") (Q\ 



-V* £ ( s in(k( Xl - 9)) - T -^ Vh {X i - 9)) 
~{ V 1 /o,fc;ll / 



This efficient central sequence can be seen as the orthogonal projection of A^C^) onto 

the subspace orthogonal to A^ fc . 1 (6') , which ensures that Ay™^^(6*) and A^ fc . 1 (0) are 
asymptotically uncorrelated. An asymptotic test can then be easily obtained by considering 
a studentized version of IS^^JJ {9) . Unfortunately, by doing so, it can be shown that the effect 
of replacing 9 with 9 is annihilated under /q only (the asymptotic linearity fails to hold for 
9o 7^ /o)- Therefore, rather than having as in the present paper a test which is valid under any 
density fo € J- with a fixed location 9 , we obtain a test which is valid under /o only (complete 
parametric test) and any value of 9 . Constructing tests which are completely distribution-free 
(with respect to both the underlying base density and the location parameter) is an ongoing 
research project. 
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A Proof of Theorem 2.1 

Our proof relies on Lemma 1 of Swensen (1985) — more precisely, on its extension in Garel and 
Hallin (1995). The sufficient conditions for ULAN in those results readily follow from standard 
arguments (hence are left to the reader), once it is shown that (9, A)' ^ (f^ x )^ 2 (x) (see (2.1)) 
is quadratic mean differentiable at any (#,0)', which we establish in the following lemma. 

Lemma A.l. Let fo&J? and k £ No, and assume that Assumption (A) holds. Define 

f k ._ 1 fo(x - 9) 



and 



D x (fl x y^(x)\ x=0 := -f^( x -9) S \n{k{x-9)). 
Then, for any 9 € [— it, tt], we have that, as (t,£) — > (0,0), 
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(i) Z, {(fe + t,o) 1/2 (x) ~ ifl^-tDetfl^ixifdx = o{t 2 ), 
(n) Z, ((^ +M ) V2 (^) " (tf +t ,o) 1/2 -lD x (f* +tx )^(x)\ x=0 ) 2 dx = o(£ 2 ), 
(m) (D x (f* +tx )^(x)\ x=0 - D x (f^/\x)\ x=0 f dx = o(l), 

h z\ui +t ^ X ) - w& -(<)'( ^;gCto ) } 2<fa = o(ll(M),||2) - 

Proof of Lemma A.l. (i) By definition of f@ we can rewrite the left-hand side of (i) under 
the simpler form 

f U% -(0 + t))- fl'\x - 9) + UJjfc^L) 2 dx. (A.9) 

Next, the a. e. -differentiability of /o (Assumption (A)) combined with the mean value theorem 
turns the latter expression into 

--t'FiP^rwrr^ 1 ,A - 10) 

J ~* V/o (a: -^*) /o (x-9)J 

with 0* € (0, + t) . Assumption (A) ensures that is continuous over the compact 

fo ( x ) 

[—ir, 7r] , hence its square can be bounded by a sufficiently large constant; consequently, the 
Lebesgue dominated convergence theorem implies that (A. 10) is o(t 2 ) . 
(ii) Similarly, the left-hand side integral in (ii) can be re-expressed as 

J f (x-e-t) (^l+£sm(k{x -9 -t))) 1/2 -1- £^sm(k(x -9 -t))^j dx. 

Exactly as for (i), the differentiability of sin(fcx) allows us to re- write this integral under the 
form 

\e £ Mx _ «, _ t) ^ {k(x _,_ t)) _ ,) 2 dx 

with £* £ (0, £). Since sin 2 (/cx)/o(x) is integrable and (1 + £* s'm(k(x — 9 — i)))" 1 is bounded 
by a constant not depending on £ (indeed, we can take £* < £ < 1/2 as £ — > 0, hence 
1 + £* sin(A;(x — 9 — t)) > 1/2 over [— 7r,7r] which does not depend on £), Lebesgue's dominated 
convergence theorem applies and yields the desired o(£ 2 ) quantity. 
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(iii) The left-hand side in (hi) equals 

\ [ [fo /2 (x-(0 + t))sm( y k(x-( y e + t)))-fQ /2 (x-e)sm(k( y x-e))ydx. (A.ll) 



1/2 

Since f (x)s'm(kx) is square-integrable, the quadratic mean continuity entails that (A.ll) 
tends to zero as t — > 0, hence is an o(l) quantity. 

(iv) The left-hand side in (iv) is bounded by C(S\ + S2 + £ 2 Ss) , where 

Si = f ((f e k +t ,o) 1/2 (x) ~ (/£o) 1/2 " tDoif^ix)) 2 dx, 
S 2 = J* ((fl +t/ ) l,2 (x) - (fl +tfi ) 1 ' 2 - £D x (f« +t>x )V 2 (x)\ x=o y dx, 

and 

53 = j\ (^A(/| +t ,A) 1/2 (x)U=0 " D x (fl x )^(x)\ X=o y dx. 

The result then follows from (i), (ii) and (iii). □ 

B Proof of Theorem 3.1. 

Proofs of Theorem 3.1. Fix /o £ T . Part (i) of the theorem trivially follows from the Central 
Limit Theorem combined with the fact that 

n *(n>,e In-^S^sinCfcpfr-fl))! \n~ 1/2 J2?=i sm(k(Xj - 9))\ 

= i -i r n ■ 27771? q\\\i/2 ~ rr U/2 Hop(l) [a.ll) 

[n E i= i sm (H x i ~ 0))) ' (r/ ,fc;22) i/2 

as n — > oo under ±u)\y.f Q - Regarding Part (ii), the situation is slightly more subtle but 
can also be readily solved by having recourse to the so-called "Third Lemma of Le Cam" 
(see Le Cam 1986). Under p|™^ i/2 («)y,^, fc , > the asymptotic normality of A^iO) with mean 

Cf (k,k')T2 and variance rf 0: k-22 is obtained by establishing the joint normality of A[?2(#) and 

log ( dP |™ ) n _ 1/2r (n )) ,. /o k ,/ dF p,oy-j ) under P (e,o)';/ and then a PP ] y in § Le Cam's third Lemma 
(which holds thanks to the ULAN property). Part (ii) follows immediately since (B.12) also 
holds under P^ , , by contiguity. Finally, Part (iii) trivially follows from (B.12) and 

(e.n-Va^ >)';f ,k' 

the optimality features of the parametric test 4>^.f for all /o £ T . □ 



References 

[1] Abe, T., Kubota, Y., Shimatani, K., Aakala, T. and Kuuluvainen, T. (2012) Circular 
distributions of fallen logs as an indicator of forest disturbance regimes. Ecol. Indie. 18, 
559-566. 



20 



[2] Abe, T. and Pewsey, A. (2011) Sine-skewed circular distributions. Stat. Papers 52, 683-707. 

[3] Azzalini, A. (1985) A class of distributions which includes the normal ones. Scand. J. 
Statist. 12, 171-178. 

[4] Azzalini, A. and Capitanio, A. (1999) Statistical applications of the multivariate skew- 
normal distributions. J. R. Stat. Soc. B 61, 579-602. 

[5] Azzalini, A. and Capitanio, A. (2003) Distributions generated by perturbation of symmetry 
with emphasis on a multivariate skew-i distribution. J. R. Stat. Soc. B 65, 367-389. 

[6] Azzalini, A. and Genton, M. G. (2008) Robust likelihood methods based on the skew-t and 
related distributions. International Statistical Review 76, 106-129. 

[7] Batschelet, E. (1965) Statistical Methods for the Analysis of Problems in Animal Orientation 
and Certain Biological Rhythms. American Institute of Biological Sciences, Washington, 
DC. 

[8] Butler, C. C. (1969) A test for symmetry using the sample distribution function. Ann. 
Math. Statist. 40, 2209-2210. 

[9] Cassart, D., Hallin, M. and Paindaveine, D. (2008) Optimal detection of Fechner- 
asymmetry. J. Statist. Plann. Inference 138, 2499-2525. 

[10] Chiogna, M. (2005) A note on the asymptotic distribution of the maximum likelihood 
estimator for the scalar skew-normal distribution. Stat. Methods Appl. 14, 331-341. 

[11] DiCiccio, T. J. and Monti, A. C. (2004) Inferential aspects of the skew exponential power 
distribution. J. Amer. Statist. Assoc. 99, 439-450. 

[12] Garel, B. and Hallin, M. (1995) Local asymptotic normality of multivariate ARMA 
processes with a linear trend. Ann. Inst. Statist. Math. 47, 551-579. 

[13] Hajek, J. and Sidak, Z. (1967) Theory of Rank Tests. Academic Press, New York. 

[14] Hallin, M. and Ley, C. (2012) Skew-symmetric distributions and Fisher information - a 
tale of two densities. Bernoulli 18, 747-763. 

[15] Hill, D. L. and Rao, P. V. (1977) Tests of symmetry based on Cramer-von Mises statistics. 
Biometrika 64, 489-494. 

[16] Jander, R. (1957) Die optische Richtungsorientierung der roten Waldameise (Formica Rufa 
L.). Z. Vgl. Physiol. 40, 162-238. 

[17] Jones, M. C. and Pewsey, A. (2012) Inverse Batschelet distributions for circular data. 
Biometrics 68, 183-193. 

[18] Jupp, P. E. and Spurr, B. D. (1983) Sobolev tests for symmetry of directional data. Ann. 
Statist. 11, 1225-1231. 



21 



[19] Kato, S. and Jones, M. C. (2012) A family of distributions on the circle with links to, and 
applications arising from, Mobius transformation. J. Amer. Statist. Assoc. 105, 249-262. 

[20] Kozubowski, T. J. and Panorska, A. K. (2004) Testing symmetry under a skew Laplace 
model. J. Statist. Plann. Inference 120, 41-63. 

[21] Le Cam, L. (1960) Locally asymptotically normal families of distributions. Univ. Calif. 
Publ. Statist. 3, 27-98. 

[22] Le Cam, L. (1986) Asymptotic methods in statistical decision theory. Springer- Verlag, New 
York. 

[23] Ley, C. and Paindaveine, D. (2009) Le Cam optimal tests for symmetry against Ferreira 
and Steel's general skewed distributions. J. Nonparam. Statist. 21, 943-967. 

[24] Ley, C. and Paindaveine, D. (2010) On the Singularity of Multivariate Skew-Symmetric 
Models. J. Multivariate Anal. 101, 1434-1444. 

[25] Ley, C, Swan, Y., Thiam, B. and Verdebout, T. (2013) Optimal R-estimation of a spherical 
location. Statist. Sinica 23, 305-332. 

[26] Mardia, K. V. (1972) Statistics of Directional Data. Academic Press, London. 

[27] Mardia, K. V. and Jupp, P. E. (2000) Directional statistics. Wiley, Chichester. 

[28] McWilliams, T. P. (1990) A distribution-free test for symmetry based on a runs statistic. 
J. Amer. Statist. Assoc. 85, 1130-1133. 

[29] Modarres, R. and Gastwirth, J. L. (1996) A modified runs test for symmetry. Statist. 
Probab. Lett. 31, 107-112. 

[30] Perez, I. A., Sanchez, M. L., Garcia, M. A. and Pardo, N. (2012) Analysis of C0 2 daily 
cycle in the low atmosphere at a rural site. Sci. Total Environ. 431, 286-292. 

[31] Pewsey, A. (2000) Problems of inference for Azzalini's skew-normal distribution. J. Appl. 
Statist. 27, 859-870. 

[32] Pewsey, A. (2004) Testing for circular reflective symmetry about a known median axis. J. 
Appl. Statist. 31, 575-585. 

[33] Rayleigh, Lord (1919) On the Problem of Random Vibrations and Random Flights in one, 
two and three Dimensions. Phil. Mag. 37, 321-346. 

[34] Rothman, E. D. and Woodrofe, M. (1972) A Cramer-von Mises type statistic for testing 
symmetry. Ann. Math. Statist. 43, 2035-2038. 

[35] Schach, S. (1969) Nonparametric symmetry tests for circular distributions. Biometrika 56, 
571-577. 



22 



[36] Swensen, A. R. (1985) The asymptotic distribution of the likelihood ratio for autoregressive 
time series with a regression trend. J. Multivariate Anal. 16, 54-70. 



[37] Umbach, D. and Jammalamadaka, S. R. (2009) Building asymmetry into circular distribu- 
tions. Stat. Probab. Lett. 79, 659-663. 

[38] van der Vaart, A. (2002) The statistical work of Lucien Le Cam. Ann. Statist. 30, 631-682. 

[39] Wang, J., Boyer, J. and Genton, M. G. (2004) A skew-symmetric representation of 
multivariate distribution. Statist. Sinica 14, 1259-1270. 



23 



